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ABSTRACT 

Consider the quadratic form f3 = y*(YY* + pl) _1 y where p is 
a positive number, where y is a random vector and Y is a TV x 
K random matrix both having independent elements with differ- 
ent variances, and where y and Y are independent. Such quadratic 
forms represent the Signal to Noise Ratio at the output of the linear 
Wiener receiver for multi dimensional signals frequently encoun- 
tered in wireless communications and in array processing. Using 
well known results of Random Matrix Theory, the quadratic form /3 
can be approximated with a known deterministic real number (3k in 
the asymptotic regime where K — > oo and K/N — > a > 0. This 
paper addresses the problem of convergence of f3. More specifically, 
it is shown here that \f~K(J3 — @k) behaves for large K like a Gaus- 
sian random variable which variance is provided. 

Index Terms — Antenna Arrays, CDMA, Central Limit Theo- 
rem, MC-CDMA, Random Matrix Theory, Wiener Filtering. 

1. INTRODUCTION 

Consider the TV dimensional received signal 
r = Ss + n 

where s = [so, si, . . . , sk] T is the transmitted complex vector sig- 
nal with size K + 1 satisfying Ess* = Ik+i, matrix S represents 
the "channel" in the wide sense and n is the independent AWGN 
with covariance matrix Enn* = pi at > 0. In this article, we are in- 
terested in the performance of the linear Wiener estimate (also called 
LMMSE for Linear Minimum Mean Squared Error estimate) of sig- 
nal so. Among the various performance indexes, we shall focus on 
the Signal to Noise Ratio (SNR) which can be expressed as follows: 
Partition the channel matrix as S — [y Yj , then the Wiener estimate 
so of so writes so = y* (SS* + pljv) _ r and the associated SNR 
(3k is given by: 

(3 K =y* (YY*+pI^)- 1 y 

A popular tool to address this problem, widely used in multidimen- 
sional signal processing and communication engineering, is repre- 
sented by the theory of Large Random Matrices: Assume that S 
is random (in this case, Pk becomes a conditional SNR) and let 
TV — > oo with K/N — > a > (denoted in the sequel by "TV — > oo" 
for short). As amply shown in the literature, there are many sta- 
tistical models related to S for which there exists a deterministic 
sequence [3k such that [3k ~ (3k — ► almost surely (a.s.); this ap- 
proximation is generally defined as the solution of an implicit equa- 
tion. Beyond the convergence of the SNR, a natural practical and 
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theoretical problem concerns the study of its fluctuations (think for 
instance to the outage probability evaluations). Despite its interest, 
there are very few related articles in the literature. In this paper, we 
provide a Central Limit Theorem (CLT) for [3k as K — > oo for a 
general model of matrix S: Assume that the TV x (K + 1) matrix 
S is given by: 



r t„ tN,k 



(1) 



where (<r^ fe ; 1 < n < TV; < k . . . , K) is a sequence of real num- 
bers called a variance profile and where the complex random vari- 
ables Wnh are independent and identically distributed (i.i.d.) with 
EW„ fc = 0, EW B \ = 0, and E|W„ fc | 2 = 1. In this case, the 
quadratic form (3k is given by: 



0K 



wjD, 1 / 2 (YY* + plj> 



K 



-l D l/2 



(2) 



where w = [Wio, W 2 q, W N0 ] 
onal nonnegative matrix Do = diag 



1 and D is the TV x TV diag- 
<Tio , • • • , u%o ) ■ An important 
special case that we shall describe carefully in the sequel is when 



the variance profile is separable, i.e a\ k = d n dk 



Among the many applications of the general model (T), let us 
mention: 



Multiple antenna transmissions with K + 1 antennas at the 
transmission side and TV antennas at the reception side. Here 
we consider the transmission model r = Hs + n where 
B = -^=HP 1/2 , matrix H is a TV x (K + 1) random 
matrix with complex Gaussian elements representing the ra- 
dio channel and P = diag(po, ■ • ■ ,Pk) is the (determinis- 
tic) matrix of the powers given to the different users. Write 
H = [ho • ■ ■ Iik], and assume the columns are indepen- 
dent, which is realistic when the transmitters are distant one 
from another. Let be the covariance matrix Cfc = Ehfch*. 
and let Cfc = UfcAfcUfc be a spectral decomposition of Cfc 
where Afc = diag((A„fc) n= i,... 1 jv) is the matrix of eigenval- 
ues. If the eigenvector matrices Uo, . . . , Uk are all equal 
(note that sometimes they are all identified with the Fourier 
TV x TV matrix |1|), then one can show that matrix S in- 
troduced above can be replaced with matrix S of Model (0 
where the W„k are standard Gaussian iid and a 2 fe = \„kPk- 
In certain situations it is furthermore assumed that Ao = 
■ • • = Ak = diag(Ai, . . . , Ajv): this is the well known Kro- 
necker model with correlations at reception. In our setting 
this model is accounted for by the separable variance profile 



= A, 



Pk- 



• CDMA transmissions on flat fading channels. Here N is the 
spreading factor, K + 1 is the number of users, and 



WP 



1/2 



(3) 



where W is the N x (K + 1) signature matrix assumed here 
to have random i.i.d. elements with mean zero, variance 1/N 
and where P = diag(po, • • • , Pk) is the users powers matrix. 
In this case, the variance profile is separable with d n = 1 and 

dk = §Pk. 

MC-CDMA transmissions on frequency selective channels. 
In the uplink, the matrix £ is written: 

£ = [Howo • • ■ Hk+iWk+i] 

where Hfc = diag(/tfc(exp(2i7r(n — l)/N)) n —i jv) is the 

radio channel matrix of user k in the discrete Fourier domain 
and W = [wi, ■ ■ • , wjc] is the N x (K + 1) signature ma- 
trix with iid elements as in the CDMA case above. Model- 
ing this time the channels transfer functions as deterministic 
functions, we have a\ h = ■^■|fefc(exp(2«7r(n — 1) /7V))| 2 . 
In the downlink, we have 



HWP 



1/2 



(4) 



where H = diag(/i(exp(2i7r(n — l)/AT)) n =i,...,jv) is 
the radio channel matrix in the discrete Fourier domain, 
the x (K + 1) signature matrix W is as above, and 
P = diag(po, ■ ■ ■ ,Pk) is the matrix of the powers given to 
the different users. Model © coincides with the separable 



variance profile case with d n 
and d k = p k . 



\h(exp(2nr(n-l)/N))\ 



About the Literature. In the communication engineering litera- 
ture, the CLT for the quadratic forms has been considered probably 
for the first time in [2 1, where the authors consider the case where £ 
is a matrix with i.i.d. elements. Their results are based on [ 3 1 where 
the asymptotic behaviour of the eigenvectors of ££* is described. 
Recently, [4| considered the more general CDMA Model l[3}. The 
model considered in this paper includes the models of |2| and |4| 
as special cases. The approach used here to establish the CLT is 
powerful yet simple. It is based on the representation of Pk as the 
sum of a martingale difference sequence and the use of the CLT for 
martingales Q. 

This paper is organized as follows. In Section[2]we recall the first 
order results that describe the limiting behaviour of 0k . The CLT for 
Pk is stated in Section[3] A sketch of proof for the CLT is presented 
in Section|4] Finally, we provide simulations in Section[5] 

2. SNR DETERMINISTIC APPROXIMATION 

Let us begin with a definition and some notations. We say that a 
complex function t(z) belongs to class S if t(z) is analytical in the 
upper half plane C+ = {z G C : im(z) > 0}, if t(z) G C+ for all 
z G C+ and if im(z)|t(z)| is bounded over C+. We introduce the 
diagonal matrices 



= diag(o- lfc ,. 
= diag(cr^ 1 ,. 



2 \ 

,°~Nk), 

2 N 



k = l, 

71 = 1, 



,N 



and the diagonal matrix functions 

T(z) = diag(ti (*),... ,**•(*)) 
T(z) = diag(ti(*), -•-,««(*)) 
that are specified by the following proposition: 
Proposition 1 ([6 7\) The system of N + K functional equations 

-1 



t„(z) = 



i k (z) = 



z {l + itr(D„T(z)) 
z (1 + £tr(D fc T(2))) ' 



1 < n < N 
I < k < K 



5{z) = itr ( D (-z(Ijv + S(z)D 
5(z) = itr ( D (-z(I K + 5(z)f> 



has a unique solution (T, T) among the diagonal matrices for which 
the t n and the tk belong to class S. Functions t n (z) and tk(z) so 
defined admit analytical continuations over C — [0, oo). 

In the separable case, we have = d^D and D„ = d n T) where 
D = diag(di, . . . , djv) and D = diag(di, . . . , dn)- In this case, 
the system described above simplifies to a system of two equations: 

Proposition 2 The system of two functional equations 



(5) 



admits a unique solution (6, 5) £ S 2 . Moreover, letting z = 
— p G (—oo,0), the system admits a unique pointwise solution 
(6(—p),S(—p))suchthat6(—p) > 0, 6(— p) > 0. 

In this particular case, the matrix functions T and T defined by 
PropositionHJare given by T = + SD)' 1 and T = -~(I + 

<5D) . The asymptotic behaviour of Pm is characterized by the 
following theorem: 

Theorem 1 f^UUHZlfJ Let Pk = — trD T(-p) where T is given 
by Proposition^ Then 



Pk~Pk 



almost surely. 



Remark 1 In matrix model ([TJ, one sometimes assumes that the 
variance profile a 2 lk is obtained from the samples of a contin- 
uous nonnegative function iv(x,y) defined on [0, l] 2 at points 
(n/N,k/(K + 1)), i.e. a\ h = n(n/N, k/(K + 1)). In this 
particular case, the sequences Pk and 5k defined in Theorem [7] 
above (and also Corollary Q] below) converge to limits that are 
solutions of integral equations (see for instance ^§][9|^. 

In the separable case, Do = doD hence we have 

Corollary 1 (^^j) Assume the separable case = d n dk- Then 

Pk x 
— dx — 

do K^oo 

where 8k = 5 with (5, 5) being the solution of System ^ at z — —p. 

Remark 2 (see also Corollary |2] below ) In the separable case, 
Pk /do often represents the SNR of user normalized to this user's 
power. Therefore, we can naturally interpret the approximation 8k 
as an asymptotic normalized SNR. This approximation, as well as 
the asymptotic variance of the normalized SNR Pk / do defined in 
Corollary\2\is the same for all users. 







3. SNR FLUCTUATIONS: THE CLT 



4. SKETCH OF PROOF 



We now come to the main result of this paper, which holds true under 
some slight technical assumptions: 

Theorem 2 Let A and A be the K x K matrices 

K 



1 itrD^D m T(-p) 2 



and 



diag 



l ,m — 1 
2 



K (l + ^trD,T(-p)) 2 
1 trD;T(-pj 

' 1=1,. ..,K j 

where T is defined by Proposition [7] Let g be the K x 1 vector 

T 



K 



trD D 1 T(-p) i 



trD DjcT(-p)' 



Then the following hold true : 

1 ) The sequence of real numbers 



Q 2 K = (E|Wio| 4 - 1)— trD^T 2 



+ ^g T (lA'-A)- 1 A- 1 g (6) 



is well defined and furthermore 



< liminf 6x < limsupO/e < oo 
K K 



2 ) The sequence (3k satisfies 



P k-Pk 
Ok 



in distribution where (3k is defined in the statement of Theo- 
rem^ 

In the separable case, one can show that Q 2 K = d^Q, K where Q. K is 
given by the following corollary: 

Corollary 2 Assume the separable case a n k = d„dk- Let 7 = 
itrD 2 T 2 and 7 = itrD 2 T 2 . The sequence 



7KE|Wio| 4 -l) + T 



p 2 77 



satisfies < liminfif Q? K < limsup K Q, K < 00, and 
in distribution. 



•AT(0,1) 



Remark 3 These results show in particular that the asymptotic 
variance Q 2 K is minimum with respect to the distribution of the W n k 
when I W n * I = 1 with probability one. In the context of CDMA and 
MC-CDMA, this will be the case when the signature matrix elements 
have their values in a PSK constellation. 



Let Q be the TV x TV matrix Q = (YY* + pl N ) \ Recall that 
the deterministic approximation of (3k is (3k = -^trDoT. Getting 
back to Equation we can write 



-j= (wSDj /2 QDj /2 wo - trDoQ 



K 

+ J= trDo(Q-T) 



K(p K -p K ) = 



act f 

= S.K + XK 

It can be shown (TO) that Eylc = 0(1/ K). On the other hand, 
by using the independence of Wo and Q and the fact that the ele- 
ments of wo are i.i.d., one can easily show that E^- = 0(1) as 
K — *• co . As a consequence, the asymptotic behaviour of \[~K((3k — 
(3k) is given by Denote by E n the conditional expectation 
E„[.] =E[.||W„,o,W„ + i,o,...,Ww,o,Y].PutE w+x [.] =E[.||Y] 

and note that Ejv+iWoDo^QDj^wo = trDoQ. With these no- 
tations at hand, we have 



zk = y,( e " 



*t\1/2/~>t\1/2 N 



<K 



The sequence Z n is readily a martingale difference sequence with 
respect to the increasing sequence of a— fields a(Y), cr(Wjv,o, Y)), 
. . . , <r(Wi,o, • • • , Wn,o, Y). The asymptotic behaviour of £k (con- 
vergence in distribution toward a Gaussian r.v. and derivation of the 
variance 9|-) can be characterized with the help of the CLT for mar- 
tingales |5 ,Ch. 35]. 

5. SIMULATIONS 

In this section, the accuracy of the Gaussian approximation is veri- 
fied by simulation. We consider an MC-CDMA transmission in the 
uplink direction. The base station detects the symbols of a given 
user in the presence of K interfering users. We assume that the dis- 
crete channel impulse response of each user consists in L = 5 iid 
Gaussian coefficients with variance 1/L. All impulse responses are 
known to the base station. 
In this case, S is given by: 

s = IvfoHoWo • ■ • y/pK+iHx+i WK4 1 ] 

where 

• Hfc = diag(/ife(exp(2i7r(n — 1)/, iV))„=i,...,jv) is the chan- 
nel matrix of user k in the frequency domain, 

• pk is the amount of power allocated to user k, 

• Wk are assumed to belong to QPSK constellation with mean 
zero and variance 1/TV. 

In this case, a 2 k is given by: 



2 



Kp k 
TV 



h k (exp(2i7r(n - 1)/TV)) | 



We denote by P the power given to the user of interest. The other 
users are arranged into 5 classes according to their powers. The 
power of each class as well as the proportion of users within this 
class are given in tabled 

Figure[T]shows the histogram of \/~K((3k~ (3k) for TV = 16 and 
TV = 64. We note that as it was predicted by our derived results, the 



Table 1. Power and proportion of each user class 
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histogram of f3 K - 3 K for N=16 histogram of /3 K for N=64 
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Fig. 1. Histogram of -jK{fi K ~ Pk) far N = 16 and N = 64. 

histogram of \/K{Pk — Pk) is similar to that of a Gaussian random 
variable. In Figure [2] the measured second moment of Pk — Pk is 
compared with Q 2 K /K. We note that convergence is reached even 
for if = 8. 

6. CONCLUSION 

The Gaussian character of the SNR at the output of the Wiener re- 
ceiver for a class of large dimensional signals described by a random 
transmission model has been established theoretically and verified 
by simulation. 
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Fig. 2. Second moment of (3k — Pk 
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